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The spin of a single electron in an electrically defined quantum dot in a 2DEG can be manipu- 
lated by moving the quantum dot adiabatically in a closed loop in the 2D plane under the influence 
of applied gate potentials. In this paper we present analytical expressions and numerical simula- 
tions for the spin-flip probabilities during the adiabatic evolution in the presence of the Rashba 
and Dresselhaus linear spin-orbit interactions. We use the Feynman disentanglement technique to 
determine the non-Abelian Berry phase and we find exact analytical expressions for three special 
cases: (i) the pure Rashba spin-orbit coupling, (ii) the pure Dresselhause linear spin-orbit coupling, 
and (iii) the mixture of the Rashba and Dresselhaus spin-orbit couplings with equal strength. For 
a mixture of the Rashba and Dresselhaus spin-orbit couplings with unequal strengths, we obtain 
simulation results by solving numerically the Riccati equation originating from the disentangling 
procedure. We find that the spin-flip probability in the presence of the mixed spin-orbit couplings 
is generally larger than those for the pure Rashba case and for the pure Dresselhaus case, and that 
the complete spin-flip takes place only when the Rashba and Dresselhaus spin-orbit couplings are 
mixed symmetrically. 



I. INTRODUCTION 

Geometric phases abound in physics and their study 
has attracted considerable attention since the seminal 
work of Berryjii^ In recent years a number of researchers 
have shown their interest in the geometric phases associ- 
ated with single- and few-spin systems for potential ap- 
plications in the field of quantum computing and non- 
charge based logic4"^ One interesting proposal is the no- 
tion that the spin of a single electron trapped in an elec- 
trostatically defined 2D quantum dot can be manipulated 
through the application of gate potentials by moving the 
center of mass of a quantum dot adiabatically in a closed 
loop and inducing a non-Abelian matrix Berry phase^S A 
recent work shows that the Berry phases can be changed 
dramatically by the applications of gate potentials and 
may be detected in an interference experiment:^ 

In the present paper, we study the non-Abelian uni- 
tary operator of the spin states during the adiabatic mo- 
tion of a single electron spin. The non-Abelian nature 
here stems from the spin-orbit coupling of an electron 
in two dimensions. The evolution operator which gives 
rise to the Berry phase is not easy to evaluate as it con- 
tains non-commuting operators. In 1951, Feynman^ de- 
veloped an operator calculus for quantum electrodynam- 
ics, in which he devised a way to disentangle the evo- 
lution operator involving non-commuting operators. In 
1958, PopovS applied the operator calculus, combined 
with group-theoretical considerations, to the spin rota- 
tion for a particle with a magnetic moment in an exter- 
nal magnetic field to obtain exact transition probabilities 
between the initial and final spin states. In a way similar 
to Popov's we employ the Feynman technique to disen- 
tangle the evolution operator for a quantum dot with 



the Rashba and Dresselhaus spin-orbit couplings and de- 
rive analytical expressions for spin transition probabili- 
ties. In particular, we obtain exact closed form expres- 
sions for three specific cases: (i) the pure Rashba spin- 
orbit couplingi^ (ii) the pure linear Dresselhaus spin- 
orbit couplingjii and (iii) the symmetric combination of 
the Rashba and Dresselhaus spin-orbit couplings. This 
approach provides us a convenient numerical scheme for 
an arbitrary mixing of the two types of spin-orbit cou- 
plings via a Riccati equation. An interesting result we 
find is that the spin-fiip probability for the case of an ar- 
bitrary mixture of the Rashba and the Dresselhaus spin- 
orbit couplings is generally greater than that for the case 
where either the Rashba or the Dresselhaus interaction 
acts alone. Furthermore, we see that the complete spin 
precession occurs only when the Rashba spin-orbit cou- 
pling and the Dresselhaus spin-orbit coupling are equal 
in strength. 

The work of Berry teaches that if parameters contained 
in the Hamiltonian of a quantal system are adiabatically 
carried around a closed loop an extra geometric phase 
(Berry phase) is induced in addition to the familiar dy- 
namical phase A slow variation of such parameters 
along a closed path C will return the system to its orig- 
inal energy eigenstate with an additional phase factor 
exp{i7„(C)}. More specifically, the state acquires phases 
after a period of the cycle T as 



|*„(T)) =exp. 



h 



En{t)dt \ ■ exp{i7„(C)} \ip„ 



(1) 

However this equation applies only to non-degenerate 
states. The detailed numerical and analytical calcula- 
tions of Berry phase 7„ (C) for the Hamiltonian of a quan- 
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turn dot in 2D plane for different non-degenerate eigen 
states are explained in Ref. [l^ The system of interest 
here (a single spin in a 2D electrically defined quantum 
dot) is degeneratei^ii^ for which ([T]) is not directly ap- 
plicable. In the formulation developed by Wilczek and 
others^'^^ for degenerate cases, the geometric phase fac- 
tor is replaced by a non-Abelian unitary operator Uab 
acting on the initial states within the subspace of degen- 
eracy. The evolution equation of the state is modified in 
the form, 

|*n,a(i)) =exp|-^^ E{t)dtjY.U-bit)\i^n,b), (2) 

where a and b are the labels for degeneracy. The non- 
Abelian unitary operator can be expressed in the form, 

Uab{t) = T exp |- J ^ Aab{t') ■ Rdt'l , (3) 
where T signifies the time-ordering, and 

Aafc = -ihilpn.a |Vr| Ipn.b), (4) 

R and Vr being a vector and the gradient in parame- 
ter space, respectively. In general, the geometric phase 
transformation Uab{t) of ([3]) in parameter space contains 
non-commuting operators and time-dependent parame- 
ters. It is possible to view the parameter-dependent evo- 
lution in the subspace of degeneracy as a non-Abelian 
local gauge transformation. Correspondingly Aab in @ 
may be seen as a non-Abelian gauge connection (or the 
Yang-Mills fields). 

Although it is not straightforward to construct the 
non-Abelian gauge connection, we consider the following 
observation instructive for the case where the parameter 
space coincides with the configuration space. Suppose 
the Hamiltonian of a system is given by 

i/=-^(P-A)2 + l/(r). (5) 
Zm 

The energy eigenequation H\ipn) — En\ipn) remains in- 
variant under the local (position-dependent) gauge trans- 
formation, 

|V;„) ^ IV^;) = C/IV-n), A^ A' = UAW +inUWUl 

(6) 

If we choose such a gauge that the transformed vector po- 
tential vanishes, that is. A' = 0, then the transformation 
operator is to be of the form, 

i7 = exp|-J j(A-dr|. (7) 

In other words, this transformation will "gauge away" 
the vector potential from the Hamiltonian ([5]). Con- 
versely, if the state with the vanishing gauge is taken 
to be the initial state, the final state with an arbitrary 
gauge A is obtained by the inverse gauge transformation, 



= U~^\iIj'^). Moreover, if the inverse gauge process 
is time-dependent via the variation of position, then the 
evolution operator is given by 

[/(t) = i7-i(t) ^Texp A -f dtj . (8) 

This observation will be useful for our discussion on the 
Berry phase associated with the spin-orbit coupling. 

A matrix element of the evolution operator gives the 
transition amplitude (propagator) from an initial state 
to the final state, which is usually evaluated by approx- 
imation. For instance, the propagator for the spin-orbit 
interaction has been calculated semiclassically in a dif- 
ferent context by Feynman's path integral represented in 
coherent stateSfi^ 

In Sec. 2, we treat the phase transformation ([3]) as a 
gauge transformation, and employ Feynman's disentan- 
gling technique, rather than Feynman's path integral, to 
evaluate the time-ordered exponential for the spin-orbit 
coupling Hamiltonian. Use of Feynman's disentangling 
method in Popov's versioiJ^ enables us to obtain analyt- 
ical and numerical results for the spin transition proba- 
bilities without approximation. In Sec. 3, we plot the 
spin-flip probability versus the rotation angle, and com- 
pare the data for the pure Rashba, the pure Dresselhaus, 
and mixed cases. Sec. 4 is devoted in deriving analytical 
expressions of the non-Abelian Berry phase (the adia- 
batic evolution operator as a 2 x 2 matrix) for the pure 
Rashba and the pure Dresselhaus coupling. 

II. SPIN TRANSITION PROBABILITIES VIA 
FEYNMAN DISENTANGLING METHOD 

To discuss the revolution of spin that induces a geo- 
metric phase, we consider a GaAs quantum dot formed 
in the plane of a two-dimensional electron gas (2DEG), 
the center of mass of which moves adiabatically along a 
closed path under the influence of applied potentialsi^ 
The single-electron Hamiltonian in 2DEG (in the xy 
plane) may be written in the form, 

H^^P^ + Hso, (9) 
2m 

where m is the effective mass. The flrst term is the kinetic 
energy in two dimensions. Evidently, = P^ + Py. The 
second term is the spin-orbit (SO) coupling Hamiltonian 
in linear approximation, 

HSO - 2aiPyS^ - P^Sy) - 2P{P,S^ - PySy). (10) 

Here S is the spin operator whose components obey the 
SU{2) algebra (see, e.g., Ref.lTgl): 

[S+,S^]^2So, [So,S±]^±S±, (11) 

where S± = Sx i iSy and 5*0 = Sz- The spin-orbit 
Hamiltonian PH)) consists of the Rashba coupling whose 
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FIG. 1. (color online) Transition probability, it)i/2,-i/2 vs. 6 for three cases: (a) pure Rashba (/? — 0), (b) pure Dresselhaus 
(a = 0), and (c) mixed (non-zero a and /?) spin-orbit interactions. The orbital radius is 60 nm. The three curves represent 
the following electric field strengths: 1 x 10^ V/cm (solid black line), 5 x 10^ V/cm (dashed red line), and 1 x 10^ V/cm 
(dotted-dashed blue line) respectively. 
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FIG. 2. (color online) Transition probability wi/2,-1/2 vs. 9 for the following cases: (a) pure Rashba (/3 — 0), (b) pure 
Dresselhaus (q = 0), and (c) mixed (non-zero a and /3). The orbit radius is chosen to be 250 nm and the following values of the 
electric field are considered: 1 x 10^ V/cm (solid black line), 5 x 10^ V/cm (dashed red line), and 1 x 10® V/cm (dotted-dashed 
blue line) 
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FIG. 3. (color online) Transition probability Wi/2,-1/2 vs. 9 for the following cases: (a) pure Rashba (/3 — 0), (b) pure 
Dresselhaus (a = 0), and (c) mixed (non-zero a and /3). The orbit radius was chosen to be 500 nm and the following values of 
the electric field were chosen: 1 x 10^ V/cm (solid black line), 5 x 10* V/cm (dashed red line), and 1 x lO*" V/cm (dotted-dashed 
blue line). 
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FIG. 4. (color online) Transition probability Wi/2,-1/2 vs. 
for the following cases: pure Rashba (/3 = 0: dotted-dashed 
blue line), pure Dresselhaus (a = 0: dashed red line) and 
mixed (non-zero a and /3: solid black line). The orbit ra- 
dius was chosen to be 250 nm and the following values of the 
electric field were chosen: (a) -E = 1 x 10^ V/cm and, (b) 
S = 5 X 10^ V/cm. 
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FIG. 5. (color online) Transition probability Wi/2,-1/2 vs. 8 
for a — j3. Physically, this situation occurs for electric field 
strength given by -E = 3.02 x 10® V/cm. Here the solid black 
line represents for both Rashba and Dresselhaus spin-orbit 
coupling effects whereas the dashed red line represents only 
for Dresselhaus spin-orbit coupling effect and open red circles 
represents only for Rashba spin-orbit coupling effect. Here we 
choose 60 nm orbit radius. 

strength is characterized by parameter a and the linear 
Dresselhaus coupling with /3. These coupling parameters 
are dependent on the electric field E of the quantum well 
confining potential (i.e., E = —dV/dz) along z-direction 
at the interface in a heteroj unction as 

e ^ ^ 0.77947c /2TOe\^^^ 

where an = 4.4 and 7c — 26 eVA^ for the GaAs quan- 
tum dotJ^ The quantum well confining potential (i.e., 
E = —dV/dz) along z-direction is not symmetric in III-V 
type semiconductor.^^ It means, the formation of quan- 
tum dot at the interface of III-V type semiconductor in 
the plane of 2DEG is asymmetric. 

Now we look for the evolution operator ([S]) for the case 
of spin-orbit coupling. It has been known that the linear 
spin-orbit term in ^ can be gauged away^^SiSi In fact. 




Rotation angle 6 (2ii) 

FIG. 6. (color online) Transition probability u)i/2,_i/2 vs. 6 
for ot = p. Physically, this situation occurs for electric field 
strength given hj E = 3.02 X 10® V/cm. The following orbit 
radii were chosen: 60 nm (solid black line), 175 nm (dashed 
red line), and, 250 nm (dotted-dashed blue line). 

the Hamiltonian ^ may be expressed as 

77 = ^ (P - A)2 - Fo, (13) 

where 
and 

Vo = mn^{a^ + P"^). (15) 

If the semiclassical momentum P = mr is used for the 
adiabatic evolution, then the spin-orbit gauge connection 
is related to the SO Hamiltonian dTU]) . 

A-Y = -Hso- (16) 

Assuming that the spin-orbit coupling is adiabatically 
introduced into the initial state, we obtain via the 
evolution operator of the form, 

C/(t) = rexp|-^^ H5o(i')rft'}, (17) 

which we shall evaluate by utilizing the Feynman disen- 
tangling method. This form of the evolution operator 
is commonly employed for Berry's phase associated with 
the spin-orbit interactioni^ii^ 

Before disentangling, we note that the SO Hamiltonian 
(flUl) may also be expressed as 

Hso=H+S++H^S^ (18) 

with 

H± = {aPy~ pP^)T^{liPy-aP^)- (19) 
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Suppose the quantum dot orbits around a closed circular 
path of radius Rq in the x — y plane under the influence 
of gate potentials, so that r = Ro{cosujt, sinwi, 0). Then 
the semiclassical momentum P = mr has components, 

Px = —Rquiuj smut, Py ~ Rquiu costjt, P^ = 0. (20) 

Substitution of ^ into (HH) yields 

H± = Ramuj{ae^"^' T j/3e=^"*). (21) 

Since 5*+ and S- do not commute, the evaluation of the 
time-ordered exponential for the evolution operator (|17|) 
is cumbersome. 

We now turn to a discussion of the Feynman disentan- 
gling technique and its application to the present prob- 
lem. For the case where the Hamiltonian is given by 

H = a{t)A + P{t)B + -f{t)C + ---, (22) 

where A, B, C, ... are noncommuting operators, and 
a, /?, 7, ... are time-dependent parameters, Feynman^ 
devised an operator calculus by which the time-ordered 
exponential can be disentangled in the form 

[/(t) - e'^W VW^e=(*)^ • • • , (23) 

where a{t), b{t), c{t), ... are time-dependent coefficients 
which can be determined by solving relevant differential 
equations. This procedure is referred to as the Feynman 
disentangling method.—™ 

Here we apply Feynman's method for disentangling the 
time-ordered exponential in (jl7p with the Hamiltonian 
(flUl) . First we rewrite the Hamiltonian ([TU)) as 

Hso = + {H+ -0S+ + , (24) 

where ^ is a time-dependent function to be determined 
appropriately. According to Feynman's procedure, the 
evolution operator may be put into the form, 

jj^^^^^a(t)s+ expjl^* dt' [{H+-OS'+ + H^S'_]Y 

(25) 

where 

a{t)^^ [\{t')dt', (26) 
in. Jo 

S'+ = 6^"'^+ 3+6"'^+ = S+ (27) 

and 

S'_ = e-^'^^ S+e'"^+ =S-- 2aSo - a^S+. (28) 

Substituting (gT]) and (US)) into (gS]) and choosing £,{t) 
such that the coefficient of in the integrand vanishes, 
we get 

U{t) = e°(*)^+ T exp 1^ ^ dt' [-2aH-So + H-S-]j , 

(29) 



in which the term containing S+ is disentangled. In a 
similar fashion, we disentangle the time-ordered expo- 
nential involving the mutually non-commuting operators 
and S- by letting 

U{t) = e°(*)^+e''(*)^° 

T exp 1^ ^ dt' [{-2aH. - 1^)3'^ + H^S'l]j , (30) 

where 

b{t) = ^ fv{t')dt', (31) 
tri Jo 

5^' = e-'''^" Soe'"^° = (32) 

and 

5// ^ e-^^'^S-.e^^" = S. e^ (33) 

Again choosing r]{t) — — 2aiJ_, we reduce the evolution 
operator (j25p into the completely disentangled form, 

j^(^) ^ga(t)S+g6(t)SogC(t)S^^ (34) 

where 

ait) = — / [H+{t') - a\t')H^{t')]dt', (35) 
in. Jo 

b{t) - - 4 f a{t')H.{t')dt' (36) 
in. Jo 

and 

1 /■* 

c{t) = — / H^{t')e^^'Ut'. (37) 
in. Jo 

Although the time-ordered exponential is disentangled, 
the evaluation of the evolution operator remains incom- 
plete until the coefficients a{t), b{t) and c{t) are deter- 
mined. In general, the integral equations ([35l) - (|37p or 
the equivalent differential equations are difficult to solve. 
In Sec. 4, we shall determine the coefficients and the 
evolution operator for the pure Rashba, and the pure 
Dresselhaus coupling. 

As it is seen in Appendix A, the spin transition prob- 
ability depends only on a{t). Therefore the full form of 
the evolution operator is not needed. To determine a{t), 
we convert the integral equation p5l) together with 
into a Riccati equation of the form, 

^^^-Ru[f[t) + r{t)a\t% (38) 

where R — mRo/h, 

f{t) = 13"^' +iae~''^\ (39) 
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and 



f*{t) = r"^* -iae'' 



(40) 



Solving psp for a(i), we can obtain the spin transition 
probabilities, Ws,s' ■ In particular, the transition proba- 
bilities from spin 1/2 to ±1/2 are calculated by 



'Wl/2,1/2 



1 



1+ a 



1+ a 



(41) 



III. NUMERICAL ANALYSIS 

As it is shown in Appendix B, exact solutions of the 
Riccati equation ([38]) can be obtained only for special 
cases, which include those for (i) the Rashba limit 
(/3 = 0), (ii) the Dresselhaus limit [a — 0) and (iii) 
the symmetric mixture of the two couplings {a — j3). 
The spin-flip probabilities obtained in Appendix B for 
exactly solvable cases (with 9 ~ ujt) are: 

(i) The Rashba limit (a 7^ 0, ^5 = 0): 



^1/2,-1/2 



^^^sin^ -yiTl]^. ; (42) 



(ii) The Dresselhaus limit (a = 0, ^5 7^ 0): 

-f/2,-1/2 = T^^«-' {^v/TTl^ q\ ; (43) 

(iii) The symmetric Rashba-Dresselhaus limit (a ~ 

/3^0): 

K/™-i/2 = |\/2ai?(sin6'- COS6I-I- 1)| . (44) 

For an arbitrarily mixed Rashba-Dresselhaus coupling 
(mixed R-D), the Riccati equation (|38|) is not exactly 
solvable. Therefore numerical analysis is needed. In the 
below we treat the mixed R-D coupling (a 7^ /3) and the 
symmetric R=D coupling (a — /3) separately. 

Comparison of the Rashba coupling, the Dres- 
selhaus coupling and the mixed R-D coupling: 

Figs. [1] [2] and [3] plot the spin-flip probability u'i/2,_i/2 
versus the rotation angle Q — ujt \a. the unit of lix for 
the orbit radius i?o=60 nm, 250 nm, and 500 nm, re- 
spectively. The plots of (a), (b) and (c) in these figures 
correspond to (a) the pure Rashba case (/? = 0), (b) the 
pure Dresselhaus case (a = 0) and (c) the mixed R-D 
case (a 7^ 0,/3 7^ 0), respectively. The three different 
values of the electric field E = \ y. 10^ V/cm, 5 x 10^ 
V/cm, and 1 x 10^ V/cm, are chosen for the curves in 
each figure, solid black, dashed red, and dotted-dashed 
blue, respectively. The symmetric case (R=D) will be 
examined separately with Figs. [5] and [6l 

The curves for (a) the pure Rashba case and (b) the 
pure Dresselhaus case are obtained from the exact results 



(l42l) and p3|) . As it is obvious from these equations, 
the spin-flip probability increases as the electric field in- 
creases via the coupling parameter but remains to be less 
than unity. Another observation we can make from these 
plots is that the periods of spin-flip for the pure Rashba 
coupling and the pure Dresselhaus coupling are different. 
This is also expected from the analytical results (j42|) and 

The curves in Figs.[TJc), [21(c) and [3l[c) show the spin- 
flip probability for (c) the mixed R-D case where both 
a and /3 are not zero and not equal. Note that they 
are not the results from the exact formula (|44| for the 
symmetric R-D coupling. Since the Riccati equation ([38]) 
for arbitrary non-zero a and /3 is not solvable, we carry 
out numerical simulations by using numerical solutions 
of ([38)) in ([4T|) . The spin- flip probability for the mixed 
case is generally larger than the pure cases. Furthermore, 
it does not reach unity if a 7^ /?. In other words, the 
complete spin-flip is not likely to occur during the entire 
period of the adiabatic motion along the closed orbit. In 
the vicinity of the symmetry point (a = /3), the transition 
probability becomes very close to unity at certain angles. 

Fig. [4] gives a further comparison study of the 
transition probability for the pure Rashba, the pure 
Dresselhaus, and the mixed case. In Fig. [D^a), when 
the electric field is weak, the curve for the mixed case 
appears to be a superposition of those for the two pure 
cases. As the electric field increases, the superposition 
effect becomes obscure as is seen in Fig. [H[b). As the 
Riccati equation is nonlinear in nature, there is no reason 
to expect that the mixed case is a superposition of the 
two pure cases. It is interesting to observe that the 
mixed case has a better chance to achieve the spin-flip 
than the pure cases during the period of evolution. 

Analysis of the symmetric R-D coupling:- The sym- 
metric mixture of the Rashba and Dresselhaus couplings 
has been discussed in connection with the persistent spin 
helixj^iS^ Bernevig et alM found an exact SU{2) sym- 
metry in the symmetric mixture and predicted the per- 
sistent spin helix which is a helical spin density wave with 
conserved amplitude and phase. Recently spin life time 
enhancement of two orders of magnitude near the sym- 
metry point {a = /3) has been reported experimentally)^ 
The coupling parameters a and /3 of the Rashba and 
Dresselhaus interactions are given by ([T^ for the GaAs 
quantum dot. The two parameters become equal at 
E — 3.02 X 10^ V/cm. For the situation in which the 
two couplings have equal strength (i.e., a = /3), the Ric- 
cati equation (j38|) is exactly solved and the corresponding 
transition probability is given by (j44p . In Fig. 5, the spin- 
flip probability versus the angle of rotation along the or- 
bit of radius 60 nm is plotted at E = 3.02 x 10^ V/cm for 
the pure Rashba case (open red circles), the pure Dres- 
selhaus case (dashed red line), and the symmetric case 
(solid black line). We see that the symmetric Rashba- 
Dresselhaus spin-orbit coupling definitely achieves a spin- 
flip during the adiabatic process whereas the two pure 
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cases have less chances. Fig. |6] plots the transition prob- 
ability of the symmetric R-D case for three different radii 
of the orbit of the quantum dot: 60 nm (solid black line), 
175 nm (dashed red line) and 250 nm (dotted-dashed blue 
line). It shows that the chance of being in the spin- flip 
state is enhanced by increasing the orbit radius. 

It is important to notice that the complete spin flip 
takes place only in the symmetric R-D coupling. This 
may be an indication of the persistent spin helix. Al- 
though the assumed orbit of motion is circular, we can 
regard the motion for a small angle of rotation as lin- 
ear. Let 6* = £ w or 6* = 37r/2 — e. If e is small, then 
sin 9 — cos 9 + 1 £j and the exact formula p4)) may be 
approximated by 



-^1/2,-1/2 



= sm 



(45) 



As e varies from to 7^/(2^/2ai^), the spin-flip probabil- 
ity moves from zero to unity, that is, the spin completes 
a full precession. For instance, if i? = 60 nm, the range 
< ^/2aRe < 7r/2 corresponds to the portion of the 
solid black curve for < 9/2tt < 0.2 in Fig. [5l Let 
Eg = tt/ (2^/2aR). Then the Rsg is the distance the elec- 
tron progresses while the spin processes by 2tt. There- 
fore, we may be able to identify this distance with the 
spin diffusion length Lg as 



Ls = Reo/n 



1 



2V2a 



(46) 



IV. 



ANALYTICAL EXPRESSION FOR THE 
NON-ABELIAN BERRY PHASE 



Applying the Feynman disentangling method, we have 
been able to reduce the time-ordered evolution opera- 
tor P7)) to the disentangled form ([M)) with the time- 
dependent scalar functions a{t), b{t) and c{t) obeying the 
integral equations (|35 |) -p7 |) . It is sometimes convenient 
to express the evolution operator as a 2 x 2 matrix in 
the spin representation of SU{2). Evidently the SU{2) 
algebra (fTTj) is satisfied by 



S+ 



1 






1 



(47) 

Using the properties S*^ = and = 1/4, we can write 
§M as 

from which immediately follows that 

This is the desired matrix representation for the Berry 
phase, and is used for calculating the spin-flip probabili- 
ties in Appendix A. 



The expressions f34|) and (|49|) remain formal until the 
time-dependent functions a(t), b(t) and c(t) are speci- 
fied. Eq. ([55]) for a{t) is equivalent to a Riccati equation 
without whose solution, and ([57]) cannot be solved 
for b{t) and c{t). In Appendix B, we show that the Ric- 
cati equation can be solved exactly if the function h{t) 
defined by 



/32 



2i?(a2 +/32)3/2 



1 



2al3 



-1 -3/2 



■ sin(2cji) 



(50) 

becomes time-independent {h{t) = ho). The last restric- 
tion ((50|l is fulfilled only when one of the following con- 
ditions is met: a = 0, /? = or a = /3. This implies that 
the function a{t) can be determined only for the pure 
Dresselhaus coupling, the pure Rashba coupling and the 
symmetric Rashba-Dresselhaus coupling. The result we 
find for a{t) is 



ait) 



if 



J4>(t) 



I/I nie*'^(*) -712^ 



where 



f{t) = /3e'"' + 



lae 



Here 



(j){t) ^ Ru{ni ~ n2) / \f{t')\dt'. 



ni,n2 = ho ± 



hi 



1. 



(51) 



(52) 



(53) 



(54) 



For convenience, we choose rii > 712. A closed form ex- 
pression for (j){t) is given in (B12). 

In calculating the spin transition probability, all we 
need is a{t). However, for completing the evolution op- 
erator we have also to determine other functions b{t) and 
c{t) by solving ([55t and (|57)) for the already determined 
function a(t). As has been mentioned above, the Riccati 
equation can be solved exactly for the pure Rashba cou- 
pling, the pure Dresselhaus coupling and the symmetric 
Rashba-Dresselhaus coupling. In the two pure couplings, 
the phase function (f)(t) can be expressed in the form. 



(55) 



where (p = + Aa^R'^ uj for the Rashba coupling and 
(p — \/l + Aji'^R'^ oj for the Dresselhaus coupling. For 
the symmetric R-D coupling, it cannot be simplified in 
the form of ((55|) . Therefore, it is difficult to carry out 
integration in ([M]) and ([57]) . This means that we have 
analytical expressions of the adiabatic evolution operator 
only for the pure Rashba and the pure Dresselhaus 
cases. For the symmetric R-D coupling, even though we 
have no analytical expressions for b{t) and c{t), we can 
calculate the spin-fiip probability since a(t) is found in 
closed form. 
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In what follows we provide the results of integration 
for the pure Rashba coupling and the pure Dresselhaus 
coupling. 

(i) The pure Rashba coupling {a ^ 0, (3^0): 
In this case, (B2), (B4), (B5) and (Ell) yield, 



2aR 



and 



(p{t) ^ (fit, (f = ujy/l + ^a^B?. 



Upon substitution of these results into (B13) we arrive 
at 



a{t) = -e" 



(56) 



where 



From Eg. (1361) . using 
together with ([56]), we obtain 



(me^^'* - r7,2)2 ^ ' 



and 



1 - e^vt 



(58) 



(ii) The pure Dresselhaus coupling (a = 0, /3 7^ 0): 
In this case, we have 



1 



and 



Hence we get 



a{t) = ie" 



(59) 



where 



Use of 



and ((59|) leads to 



{ni - n2)2e*(^+")* 
(me*'^* - 712)2 ' 



(60) 



and 



cit) 



iipt 



nie 



n2 



(61) 



V. CONCLUSION 

In the present paper we have considered spin manip- 
ulation via the non-Abelian Berry phase induced by 
an adiabatic transport of a single spin along a circular 
path in the 2D plane in the presence of the Rashba 
and Dresselhaus spin-orbit couplings. We have adopted 
the Feynman disentangling technique to calculate the 
spin-flip probability. We have shown that the problem 
can be solved exactly in three cases: (i) the pure Rashba 
coupling, (ii) the pure Dresselhaus coupling, and (iii) 
the symmetric combination of Rashba and Dresselhaus 
couplings. For an arbitrary combination of the two 
couplings, we have carried out numerical simulations. 
We have plotted the spin-flip probability versus the 
angle of the adiabatic rotation with various values of 
the electric field and the radius of the circular path 
in the 2D plane. We have observed that a complete 
spin-flip (a complete spin precession) occurs only when 
the strength of the two couplings becomes equal. The 
relation between the complete spin precession and the 
persistent spin helix will be discussed in detail elsewhere. 
We have also obtained analytical expressions of the 
non-Abelian Berry phase for the pure Rashba case and 
the pure Dresselhaus case. 



Appendix A: The spin transition probabilities 

Following Popov's procedure, we show that the 
spin- flip probability can be expressed in the form of (|41l) . 
Since the time-evolution of the spin state can be achieved 
by a time-dependent rotation, the transition amplitude 
for spin a to a' is given by 

(a|C/(i)|a') =Z?^.,,(^,t?,0) =exp[-^(a^ + a»]d,,.(^^)• 
(Al) 

Here ip, <j} are the time-dependent Eulerian angles, 
r*^ g./ {tf, (f) are the elements of the Wigner Z?-matrix 
being the irreducible unitary representations of SU (2) 
group, and do-o-' (^9) is Wigner's d-function. 

The corresponding transition probability along the z- 
axis is 



Wo 



\d%Am)? 



(A2) 
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In particular, the transition probability from spin 1/2 to 
±1/2 is 



and 



because 



,1/2 



Wl/2,-1/2 



V 2 ; ' 



(A3) 



(A4) 



V2,i/2W =^°^§' ^^V2,-i/2W =«sin-|. (A5) 

For spin s — 1/2, the rotation matrix is given in the 
standard formJ^ii^ 



D(^,z?,V) 



where 



a = cos — exp 
2 



.ijj + (fi 



/? a* 



, /3 — i sin — exp 



(A6) 



.tp - if' 



Comparison of the evolution operator for the spin 1/2 
transition expressed in the matrix form, 



U 



and the rotation matrix yields 



ar = tan^ 



(A8) 



(A9) 



Again comparing this result with (A3) and (A4), we ar- 
rive at 



Note that wi/2,1/2 + ^1/2,-1/2 = 1- 



(AlO) 



This equation contains the Rashba limit (a ^ 0, /3 = 0), 
and the Dresselhaus limit (a = 0, /3 ^ 0), both of which 
have exact solutions. 

First we let a{t) = g{t)X{t) in (Bl). If we further let 
g{t) = — j//|/|, then we see that X{t) obeys 



'^=iRw\f{t)\{X^ ~2h{t)X^l] 



where 

and 

h{t) 



|/(t)| = [a^ + + 2a/3sin(2wt)] 



1/2 



(B3) 



(B4) 



2i?(a2 + /32)3/2 



2a[i 
q;2 + 



-1 -3/2 



■ sin(2wt) 



(B5) 

Now we consider a special case where h{t) is a constant, 
say, ft-o- In this case, (B3) can be expressed as 



dX 



(X-ni)(X-n2) 
where ni and n2 are roots of 

X"^ ~ 2hoX -1 = 0, 

that is. 



iRuj\f{t)\dt, (B6) 



ni, n2 



^ho±\ hl + 1. 



(B7) 



(B8) 



Note that 



ni«2 = -1, ni+n2^2ha, ni - 712 = 2y /ig + 1. 

(B9) 

Upon integration, we obtain with the condition X{0) = 0, 



Xit) 



I _ g«0(t) 



712 — 'T'le 



(BIO) 



The phase function (j){t) is 



Appendix B: Special solutions of the a — j3 Riccati 
equation 

Here we wish to solve under a special condition the 
Riccati equation (p8)) : 



da{t) 
dt 



^-Ruj{f{t) + r{t)a'it)}, (Bl) 



where 



(B2) 



cl>{t) = Ruj{ni - n^) / |/(r)|dr (Bll) 



which can be expressed in closed form. 



<l){t) = 2Ruj^hl + l(a + /3) 

TT 2^/a;5 
4' a + 



TT 2^/aP 
'4' a + /3 



,(B12) 



where -E(<p, fc) is the elliptic function of the second kind 
defined by 

E{ip,k)^ r Vl - fc2 sin"^ 9 d9. 
Jo 
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Consequently, for the case where h{t) = ho, the start- 
ing Riccati equation (Bl) is exactly solved, the result 
being of the form, 



a(t) = 



if e'"^ - 1 
|7|nie'«^-n2' 



(B13) 



Since 0(0) = 0, it is evident that a(0) = 0. Using (B8) 
in (B13), we obtain 



hl + 1- sin^((^/2) 



(B14) 



The transition probabilities from spin 1/2 to ±1/2 are 
given by 



W'l/2,1/2 



1 



1 



/i^ + l V 2 



and 



hi 



. 2 f<f'it) 



(B16) 



which are characterized only by the constant ho and the 
phase function (j){t). 

Although the above results are exact under the as- 
sumption that h{t) = /lo is a constant, they are approxi- 
mate results when h{t) w ho- 

Finally, specifying the values of ho and 4>{t), we shall 
obtain the exact results for the Rashba, the Dresselhaus 
and the symmetric cases. 

(i) The Rashba limit (a ^ 0, /3 = 0): In this case, 
from (B5) follows 



hQ = - 



2aRuj 



(B17) 



Furthermore the right-hand side of (Bll) can be easily 
integrated, so that 



(j,{t) = ^/l + Aa^R'^ujt. (B18) 
Thus the spin flip probability is obtained in the form. 



w 



R 

1/2,-1/2 



(B19) 



where 9 = wt. 

(ii) The Dresselhaus limit (a = 0, /3 7^ 0): In this 
case, (B5) leads to 



ho 



1 



(B20) 



The integral of (Bll) yields 



Hence the spin-flip probability is 



(B21) 



(B15) u;f/2,_i/2 - Y^^l^ sin2|iv/l + 4/32ii2^| . (B22) 



(iii) The symmetric case (a = /3 ^ 0): In this partic- 
ular case. 



ho = 0. 



(B23) 



The phase factor becomes 



0(t) = 2^/2aR [sm{ojt) - cos(wi) + 1] . (B24) 



The corresponding spin-flip probability as a function of 
6 = ojtis 



w 



1/2,-1/2 



= sin2|\/2ai?[sin6'-cos6'-M]|. (B25) 
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